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ABSTRACT 

Context. We present the basic ideas and first results from the project we are carrying out at present, the search for and characterisation 
of pre-main sequence (PMS) stars among the members of Galactic young clusters. The observations of 10 southern clusters, nine of 
them located in the Carina-Sagittarius spiral arm of the Milky Way are presented. 

Aims. We aim at listing candidate PMS member stars in young clusters. The catalogued stars will serve as a basis for future spectro- 
scopic studies of individual objects to determine the properties of stellar formation in the last phases before the main sequence stage. 
Properties such as the presence of residual envelopes or disks, age spread among PMS members, and the possible presence of several 
episodes of star formation in the clusters, are to be addressed. 

Methods. Multicolour photometry in the UBVRcIc system has been obtained for 10 southern young clusters in the fourth Galactic 
quadrant, located between Galactic longitudes 1=238° and 1=310°. For six clusters in the sample, the observations presented here 
provide the first published study based on CCD photometry. A quantitative comparison is performed with post-MS isochrones, and 
PMS isochrones from three different evolutionary models are used in the photometric membership analysis for possible PMS stars. 
Results. The observations produce photometric indices in the Johnson-Cousins photometric systems for a total of 26962 stars. The 
matching of our pixel coordinates with corresponding fields in the 2MASS data base provides astrometric calibration for all cataloged 
stars and JHK 2MASS photometric indices for 60 % of them. Post-MS cluster ages range from 4 to 60 Myr, whereas the photometric 
membership analysis assigns PMS membership assignment to a total of 842 stars, covering an age range between 1 and 10 Myr. This 
information on the PMS candidate members has been collected into a catalogue, named DAY-I, which contains 16 entries for 842 
stars in the field of 10 southern clusters. 

Key words, open clusters and associations - stars: pre-main sequence - astronomical data bases: catalogues 



^ ■ 1. Introduction These observations are feasible with small telescopes and can 

be obtained for a wide sample of clusters in the nearest corn- 
Stars clusters are commonly accepted as the astrophysical ob- plexes of star f orma tion, such as Carina and Cygnus, as well as 
jects best-suited to achieving reliable ages and distances of mem- for some c i usters i ocat ed at greater distances in the direction of 
ber stars at any location in the Milky Way. Furthermore, compar- me G a j ac tic anticentre 

ing photometric colour-magnitude (CM) diagrams with model These obser vations typically refer to older stars than those 
isochrones simultaneously provides a way to estimate masses located in embedded regions of star formation, such as the Orion 
and chemical compositions and to set constraints on the physics nebula . But they stm provide invaluable information on the ear- 
of the models themselves. In particular, the systematic observa- Uest phases of star formation . In particu i ar , the y are expected to 
tion of the Galactic sample of young clusters offers a number of contain clues to me mass distribution of forming stars in selected 
interesting possibilities for improving knowledge of how the star regionS; and they can a i so be studied with advantage at optical 
formation occurs on large scales in the Galaxy (e.g. Alfaro et al. wave i e ngths. In this spectral region, the measurement of distance 
1991) and how stars form inside the clusters (e.g. Lada 2005, and age by means f ZAMS and post-MS isochrone fittings to 
Lada & Alves 2004). By young clusters we understand objects the brightest part of the colour-magnitude diagram (CMD) is 
with ages between 1 and 10 Myr, which can be easily and pro- more reliable and leads tQ mor£ accurate estimates of physical 
ductively observed at both optical and infrared wavelengths and paramet e r s for the member stars than for isolated objects, 
with telescopes around 1 m in diameter. For objects located not pMS stars fa cmsters betwe en 1 and 10 million yr are es- 
further than 3-4 kpc, we are left with clusters containing observ- m interesting for investigating whether O-star formation 
able PMS members, also called Post-TTaun stars (Mamajek et terminate s or , on the contrary, triggers the formation of less mas- 
al. 2002), with spectral types from A to K, detectable in photo- siye ^ afj weU as the possible defidt - m the latter case( s un g et 
metric diagrams deep down to V 21-22, depending on reddening. al 1998) Furthermor e, they trace the evolution through the last 
phases of circumstellar disk dissipation and angular momentum 
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2001; Slesnick et al. 2002; Sicilia-Aguilar et al. 2005; Hartmann 
2005). In particular, the proportion of CTTs to WTTs, possi- 
bly as a function of age, can be investigated (Herbig 1998). 
Recent work reveals the existence of long-lived disks around 
low-mass PMS stars, favourable for the formation of planetary 
disks (Lawson et al. 2004; Lyo et al. 2004), and a decrease in the 
time scale of disk dissipation with increasing mass (Hernandez 
et al. 2005). 

The presence of PMS members is addressed in most stud- 
ies devoted to one or more young clusters, but specific mem- 
bership assignments are present only for close or very bright 
objects. Ever since the first systematic investigation by Walker 
(1961, and references therein) of a group of bright young galac- 
tic clusters (NGC 2264, NGC 6530, IC5146, NGC6611), the 
presence of PMS stars among the members of young clusters 
has been mentioned frequently in the literature. Many of these 
studies mention only the location of some stars on photometric 
diagrams, suggesting their PMS nature (FitzGerald et al. 1978; 
Forbes & DuPuy 1978; Sagar & Joshi 1979, 1981; Stone 1988; 
Forte & Orsatti 1984; Witt et al. 1984; Heske & Wendker 1984; 
Verschueren et al. 1990; Marschall et al. 1990; Munari et al. 
1998; Carraro et al. 2001; Patat & Carraro 2001). 

A few studies include comparison with models and eventu- 
ally more detailed observations of the most conspicuous PMS 
candidate members. Their main interest is focused on indica- 
tions of coevality or non-coevality of the cluster members, the 
properties of the cluster mass function, and the possible pe- 
culiarities of the spatial distribution of the PMS member stars 
(Stauffer et al. 1989; Trullols & Jordi 1997; Marco et al. 2001; 
DeGioia-Eastwood et al. 2001; Marco & Negueruela 2002). In 
some cases, membership probability is assigned on the basis of 
statistical considerations or discussed in some detail based on 
the global appearance of the CM diagrams as compared to PMS 
isochrones (Baade 1983; Baume et al. 1999; Baume et al. 2003; 
Prisinzano et al. 2005) 

On the other hand, detailed studies of the PMS candidate 
members, included in extensive analyses of the parent clusters, 
have been published mainly for a reduced collection of clusters 
(The et al. 1985; Hillenbrand et al. 1993; de Winter et al. 1997; 
van den Ancker et al. 1997; Sung et al. 1997; Herbig 1998; Sung 
et al. 1998; Baume et al. 1999; Delgado et al. 1999; Sung et al. 
2000, van den Ancker et al. 2000; Park et al. 2000; Park & Sung 
2002; Delgado et al. 2004, 2006). Some of these studies use pre- 
viously established PMS members, while others make their own 
PMS membership assignments, whereby the methods include lo- 
cation in several photometric diagrams, considerations of red- 
dening, spectroscopic information for some stars, and values of 
specific photometric indices such as (R-Ha). 

Most recently, several published works contain exhaustive 
analyses of the PMS populations in regions close to the Sun, 
based on the widespread data sets available and including new 
observations (Mamajek et al. 2002; Luhman et al. 2003; Sartori 
et al. 2003; Pozzo et al. 2003; Lyo et al. 2004; Prisinzano et 
al 2005). For these intensively studied objects, matches with 
available proper motions catalogues and lists of associated X- 
ray sources are used, and they provide the best-studied samples 
of PMS stars. Furthermore, in some special cases - clusters lo- 
cated at close distances, although not necessarily younger than 
10 Myr- detailed observations have been recently carried out to 
compare ages determined from isochrone fits, both for post-MS 
and for PMS models, with ages derived from the theoretical pre- 
dictions on lithium depletion during the PMS evolution (Jeffries 
& Oliveira 2005; Palla et al. 2005). 



The conclusions about age spread, mass functions, and dif- 
ferences between evolutionary models differ between studies. 
There is a general trend to admit the existence of a significant 
age spread, which also produces age differences as estimated 
from PMS and post-MS isochrones. Values suggested for age 
spread among PMS candidate members are roughly 5 Myr (Sung 
et al 2000; Park et al. 2000; Park & Sung 2002), although val- 
ues from 1 Myr (Hillenbrand et al. 1993) to more than 10 Myr 
(DeGioia-Eastwood et al. 2001; Pozzo et al. 2003; Prisinzano 
et al. 2005) are claimed. On the other hand, results from dif- 
ferent analyses of the closest associations, using a great wealth 
of available data offer different conclusions about the age differ- 
ences between association members; see Mamajek et al. (2002) 
and Sartori et al. (2003) on the age spread in the Sco-Cen OB as- 
sociation. With respect to the mass function, some indications of 
a dependence of the results on the physical conditions of star for- 
mation is found in recent analyses (Luhman et al. 2003; Pozzo 
et al. 2003), whereas recent observations of PMS stars in the in- 
termediate mass range in 30 Dor (LMC) suggest no important 
differences between star formation in clusters in the Milky Way 
and outside it (Brandner et al. 2001). 

In this paper we present the results for ten southern clus- 
ters, mostly located in the Carina-Sagittarius spiral arm of the 
Milky Way. This research is the continuation of a project, cen- 
tred on the issues of interest already mentioned in this section, 
which has produced results for some clusters, mainly located in 
the Cygnus region and surroundings (see Delgado et al. 2004, 
and references therein). The main point in this research consists 
in assigning membership to individual candidate PMS stars on 
the basis of theoretical isochrone fitting on various CMD, in the 
same way as classical multicolour photometry is used to estimate 
membership of main sequence (MS) stars by means of ZAMS 
fitting. The procedure has been shown to perform better, when 
the reddening of the parent cluster is lower (Delgado & Alfaro 
2000). The project aims at obtaining homogeneous sample of 
PMS stars at different locations in the Galactic disk. It should 
provide a sample for future spectroscopic studies of individual 
objects, as well as for analysing (i) key issues of star formation 
in clusters, such as age and relative spatial distributions of MS 
and PMS members and (ii) the presence and dissipation time 
scales of circumstellar disks. Other points of interest are the as- 
sessment of an adequate T^^-colour calibration for PMS stars 
and determining to what extent the application of calibrations 
valid for dwarfs is adequate for these stars (Delgado et al. 1998; 
Mamajek et al. 2002). Very interestingly, Lyo et al. (2004) pro- 
vide an observational isochrone from 77-Cham observations, al- 
though for PMS stars with spectral types later than K7. 

Section 2 contains a description of the observations and of 
the reduction and calibration procedures. A comparison with 
the previously published photometry of some of the clusters is 
presented. In Sect. 3 we describe and discuss the methodology 
used for determining the main physical properties of the host 
clusters, as well as the membership assignments and ages for 
post-MS and PMS stars in the observed area. Matches with the 
Two Micron All Sky Survey (2MASSf] catalogue and the result- 
ing astrometry and near-infrared colour indices of the candidate 
members are addressed. Section 4 contains a brief summary of 
how the catalogue has been elaborated and a description of its 
contents. 
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2. Observations 

The observations were carried out with the YALC0 tele- 
scope, operating at Cerro Tololo Observatory (CTIO), with the 
ANDICAIVQ camera, in two campaigns during December 2000 
and May 2001. The optical channel of ANDICAM incorporates 
an array of 2024x2024 pixels, with a plate scale of 0.3 "/pixel. 
The BVRcIc filters were mounted, together with the Stromgren 
ultraviolet u filter, instead of the commonly-used U passband 
from the Johnson system. 

Several frames were secured through each band for every 
cluster field (Table 1), with integration times for short exposures 
of 5 s for BVRI bands and 10 s for the u band, and long expo- 
sures of 1200, 600, 300, 300, and 700 s, for the uBVRI bands, 
respectively. Three fields from the catalogue of standard stars 
by Landolt (1992) were also observed. The fields of SA98-650, 
and Rubin 149, observed in the 2000 campaign, include 24 stan- 
dard stars. The fields from the 2001 campaign, SA1 10-229, and 
PG1657+078, contain 8 standards. All these fields were ob- 
served every night at a minimum of three separated airmasses, in 
every band. As mentioned above, the ultraviolet filter at our dis- 
posal was the u Stromgren filter. No measurements of the stan- 
dard fields were secured through this filter. The calibration of the 
(U - B) colour index is feasible if we dispose of accurate pub- 
lished Johnson (U - B) values in the field, and a reliable zero- 
point correction between instrumental and "standard" values can 
be determined. 

The different frames were reduced with the routines in 
the IRApQ package. Aperture corrections and atmospheric ex- 
tinction corrections were applied to the magnitudes in each 
frame. Average extinction coefficients for CTIC0 were used. 
Instrumental magnitudes were then calculated as weighted av- 
erages of the magnitudes in all frames. The weights used are the 
inverse of the squared PSF photometric errors. 

The standard calibration of UBV indices was carried out by 
direct correlation of instrumental values with published pho- 
tometry for the individual clusters, included in the WEBDA0 
database. These correlations produce accurate values for our pur- 
pose and offer a most reliable way of obtaining well-behaved 
UBV indices (Delgado & Alfaro 2000). The calibrations consist 
in determining zero-point corrections from the comparisons be- 
tween instrumental and published values, with eventual colour- 
dependent terms after examining the residuals. The root-mean- 
squared deviation of the final residuals is adopted as the uncer- 
tainty on the calibration. The calibration of the (V - R)c and 
(V - I)c indices was performed with the help of the observed 
Landolt standard stars, using the classical linear model for the 
calibration equations (see Delgado et al. 1998). As before, the 
root-mean-squared deviations of the fits were taken as the un- 
certainty on the calibrations. 

The total errors assumed for the final standard indices are 
quadratic sums of the photometric errors and calibration uncer- 
tainties as described above. In Table 1, we list the target clusters' 

2 YALO is the Yale-AURA-Lisbon-Ohio consortium (Bailyn et al. 
1999). http://www.astronomy.ohio-state.edu/YALO/ 

3 http://www-astronomy.mps.ohio-state.edu/ 
~depoy/research/instrumentation/andicam/andicam.html 

4 The Image Reduction and Analysis Facility (IRAF) is distributed 
by the national Optical Astronomical Observatory, which is oper- 
ated by the Association of Universities for Research in Astronomy, 
Inc. (AURA) under cooperative agreement with the National Science 
Foundation 

5 http://www.ctio.noao.edU/facil/s4/sec4.html#b 

6 http://www.univie.ac. at/webda//webda.html 
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Table 1. Observed clusters 



Cluster 


RA 


DEC 




1 


b 


Ref. 1 




Epoch 


2000 






dc<. 


r 




NGC 2367 


7 20 04.54 


-21 


53 


02.7 


235.65 


-3.84 


a 


NGC 3293 


10 35 48.77 


-58 


13 


28.1 


285.85 


0.07 


b 


Collinder 228 


10 43 01.26 


-60 


00 


44.8 


287.52 


-1.03 


c 


HogglO 


11 10 55.02 


-60 


22 


18.1 


290.80 


0.10 


d 


Hogg 11 


1111 32.63 


-60 


22 


49.3 


290.89 


0.14 


c 


Trumpler 18 


1111 26.29 


-60 


40 


19.1 


290.99 


-0.14 


f 


NGC 3590 


11 12 55.22 


-60 


46 


30.7 


291.21 


-0.18 


d 


NGC 4103 


12 06 42.65 


-61 


14 


26.8 


297.57 


1.18 


g 


NGC 4463 


12 29 57.47 


-64 


47 


55.7 


300.65 


-2.01 


h 


NGC 5606 


14 27 47.30 


-59 


38 


34.8 


314.84 


0.99 


i 



1 a: Vogt & Moffat 1972. b: Turner et al. 1980. c: Carraro & 
Patat 2001. d: Moffat & Vogt 1975. e: Claria 1976. f: Vazquez 
& Feinstein 1990. g: Wesselink 1969. h: Moffat & Vogt 1973. i: 
Vazquez & Feinstein 1991. 



Table 2. Comparison to published CCD photometry 



Cluster (V) (U - B) (B - V) (V - R) {V - I) Ref. 1 



NGC 3293 


0.00 


-0.04 


-0.01 


-0.01 


-0.00 


a 




19 


11 


11 


21 


10 




NGC 4103 


0.04 


-0.03 


-0.02 


0.04 


-0.03 


b 




13 


13 


13 


12 


14 




NGC 4103 


0.04 
9 




0.06 
12 




0.04 
8 


c 


NGC 5606 


0.02 


0.01 


-0.02 


0.06 


0.05 


d 




12 


11 


16 


10 


9 




Collinder 228 


-0.02 
11 


-0.01 

8 


0.00 
9 


0.03 
8 


0.04 
11 


c 



1 a: Baume et al. 2003. b: Sagar & Cannon 1997. c: Sanner et al. 
2001.d: Vazquez et al. 1994. e: Carraro & Patat 2001. 



equatorial and galactic coordinates and the photometric study of 
reference used to compute the UBV calibration. In Fig. 1 we rep- 
resent the variation in the uncertainties in a joint plot with V for 
all the measured stars together. The different levels, appreciable 
at the brightest magnitudes, reproduce the different uncertainties 
of the calibration for different clusters. 

A comparison between our absolute photometry and those 
published for the clusters with CCD data, provides the values 
listed in Table 2. In this table we list the average differences 
in the first row, and the root-mean-squared (rms) deviations for 
each cluster in hundreths of magnitude in the second row. The 
high dispersion in some cases originates in the large errors for 
faint stars. We do, however, note the absence of significant sys- 
tematic differences between our adopted colours and the pub- 
lished ones. Restricting the comparison to stars brighter than 
V=20, the rms deviation of the differences decreases below 0.1 
in all cases. 

For Collinder 228, the U BV calibration in terms of published 
photoelectric values (Feinstein et al. 1976; Turner & Moffat 
1980) provides unreliable results. Just for this cluster, the CCD 
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Fig. 1. Variation with V in the total uncertainties of standard indices for all observed clusters 



values by Carraro & Patat (2001) were then used to calibrate the 
V, (U - B), and (B - V) indices. 



3. Construction of the catalogue 

Our main objective is to detect "bona fide" PMS cluster-member 
candidates from multicolour photometry in the area of the clus- 
ter. Candidates not only need to be assigned as cluster members, 
but also have to be located in a region of the CM diagrams, which 
are compatible with a PMS evolutionary stage. We therefore de- 
vote some room to explaining the methodology used to assign 
cluster membership to the individual stars and to distinguishing 
between post-MS and PMS members. This information can be 
obtained from the comparative analysis of the observational CC 
and CM diagrams, and the comparison between observational 
and model diagrams. 



3.1. Colour excess and distance: MS membership 

The first step in the procedure of membership assignement is 
the calculation of color excess and distance for MS stars. These 
values are to be used as reference to estimate the membership 
for PMS candidate-members. We note that most of the clusters 
in our sample do not have membership assignments for stars 
in the the region of the CMD where PMS members would be 
expected. Only NGC 3292, the most-studied object, has mem- 



bership assignments for stars fainter than My =2 (Baume et al. 
2003). Some of our observed stars in Collinder 228 could be in- 
cluded in the BVRIHa photometric study of a stellar field in the 
Carina region by Hagele et al. (2003). However, this study does 
not provide published or accessible data that could be compared 
to our results. 

To calculate colour excess and absolute magnitude from the 
(U - B), (B - V) CC diagram, we used the standard values of 
the reddening parameters, a = E(U - B)/E(B - V) = 0.72, 
R = A V /E(B - V) = 3.1. For NGC 3293, the detailed study of 
the reddening slope in the cluster field by Turner et al. (1980), 
suggests the value a = 0.76, so it has been adopted here. 



3.1.1. Selection of probable, unevolved MS members 

Our procedure for first selecting the upper MS probable mem- 
bers is to use the location of stars in the CM and CC diagrams, as 
described in Delgado et al. (1998). We start with the assumption 
that all clusters have a certain number of unevolved MS mem- 
bers with spectral types bluer than AO. A tentative selection of 
these stars is made visually in the V, (B - V) CM diagram, also 
checking that the stars are congruently located in the other CM 
diagrams. This procedure is shown in Fig. 2, where the example 
of NGC 2367 is plotted. For the stars tentatively selected (left 
panel) a colour excess and a corresponding absolute magnitude 
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are calculated by means of the calibrated ZAMS line (Schmidt- 
Kaler 1982). 

The quantity V - 3.1 x E(B - V) is then plotted versus My, 
also considering previously published values of cluster param- 
eters. In the right panel of Fig. 2, we see this plot for the ten- 
tative member stars, with the two different published values for 
the distance modulus of the cluster indicated by two different 
lines (Vogt & Moffat 1972; Kharchenko 2005, also adopted in 
WEBDA). In view of this plot, we select stars as probable un- 
evolved MS members and others as evolved post-MS members, 
which will be used later to perform a quantitative estimate of the 
cluster age. The selection performed in these plots is based on 
how stars that evolved away from the ZAMS will show a lower 
distance modulus than the unevolved ones. In the same way, ro- 
tating or binary stars would have the same effect upon the DM 
calculation so are also rejected, as for the star marked with a 
triangle in the right panel of Fig. 2. 

For most clusters in our sample, this procedure provides a re- 
liable selection of upper non evolved MS members and post-MS 
members. For some of them, the selection of a cluster sequence 
is not unambiguous. In Fig. 3 we show the case of Trumpler 18 as 
an example, where two associations seem to overlap in the same 
field, at different distances. The selected possibility, according 
to previously published values, is shown with the larger distance 
modulus that would account for a star sequence located further 
away. 

The selection of probable unevolved MS members is obvi- 
ously of central importance for the membership analysis. The 
solutions adopted here, based on the plots like the one in Figs. 2 
and 3 are aimed at obtaining a basis for further PMS member- 
ship analysis, and they follow the appearance of our own CM di- 
agrams, together with previously published results. On the other 
hand, there are interesting features in the CM diagrams of several 
clusters, suggesting the presence of several groupings at differ- 
ent distances in the same line of sight. This fact is apparent for 
lines of sight crossing the Carina spiral arm at different positions. 
These features can be guessed, for example, in the plots shown 
in both Figs. 2 and 3. In this paper we only elaborate a catalogue 
of candidate PMS cluster members. The discussion of these fea- 
tures in the CM diagrams of the target clusters is postponed for 
a future paper. 

The method is applied to every cluster, and provides aver- 
age values of colour excess and distance modulus for the prob- 
able unevolved MS members. These are listed in Table 3 with 
their rms deviations, together with the results derived by other 
authors. The values in Table 3 agree in most cases with the pub- 
lished results. There are, however, some appreciable differences, 
as in the case of Hogg 10, where the three published values of 
colour excess are noticeably higher than our value. This differ- 
ence is due to our restricting the selection of probable unevolved 
MS members to avoid including binaries or evolved stars as 
much as possible. But we also generally obtain a more complete 
coverage of the cluster MS, whereas the photoelectric studies 
used contain only a few stars. In the particular case of Hogg 10, 
for example, these are the ones with highest reddening. 

3.1 .2. Membership assignment for possible MS members 

The average values of colour excess and distance modulus for 
the probable unevolved MS members are used to estimate mem- 
bership for the remaining stars. 

Values of colour excess and distance modulus are calculated 
for each star by comparing their indices with the ZAMS ref- 
erence line (Schmidt-Kaler 1982) in the (U - B),(B - V) CC 



diagram. The representative point of each star is shifted in the 
CC diagram, with the adopted reddening slope, and surrounded 
by an ellipse of semiaxes given by the uncertainties in (U - B) 
and (B - V) (Fig. 1). The minimum distance at every intersection 
of this ellipse with the reference line provides a coluor excess 
E(B - V) and a particular value of the absolute visual magnitude 
My, wich are used to calculate a value for the intrinsic distance 
modulus. 

The calculation is visualised in Fig. 4, using three stars in the 
cluster NGC 3590 as an example. Figure 4 shows how more than 
one solution is possible for each star, depending on its location in 
the CC diagram. A star is then considered as a member if any one 
of the pairs of values E(B - V), DM coincides, within the errors, 
with the average values from the probable unevolved MS mem- 
bers. In this comparison, the errors used are squared quadratic 
sums of the uncertainties plotted in Fig. 1, and the rms deviation 
of the average values (Table 3). In practice, when selecting MS 
members, we require that the colour excess be larger than the 
average value of the probable unevolved MS members minus 
the error and smaller than this value plus twice the error. When 
selecting PMS members, we simply require the color excess be 
above the value of the probable unevolved MS members minus 
the error. We try in this way to favour the detection of possible 
members, PMS members in particular, which are more affected 
by reddening than the probable unevolved MS members. 

We note that a star may already have one calculated value for 
the colour excess that fulfils the requirements for membership, 
prior to the distance modulus comparison. In considerations be- 
low we refer to these stars as reddening candidates. In fact, this 
constitutes the first membership filter. Obviously, all stars as- 
signed as members are first reddening candidates, but only stars 
that fulfil this condition and the agreement of the corresponding 
distance modulus are kept as candidate members. 

3.1 .3. Age estimate for post-MS stars 

The age of the cluster is estimated from comparing of the CM di- 
agrams with theoretical post-MS isochrones. A simple quantita- 
tive procedure is used, which provides a systematic age estimate 
instead of the commonly-used guess from a visual comparison 
between data and isochrones. A similar procedure has been used 
byTrullols&Jordi(1997). 

As done for the probable unevolved MS members, the post- 
MS probable members used in the calculations are selected from 
plots like those shown in Figs. 1 and 2, also considering the con- 
sistent location of every star in all CM diagrams, and helped by 
results from previous studies. As is known, in addition to ob- 
servational errors and to metallicity variations, which are not ex- 
pected in the studied young clusters, the presence of rotation and 
binarity cause a broadening of the cluster's main sequence (see 
for instance Fernandes et al. 1996). This would lead to overesti- 
mating the cluster age. In our selection we try to avoid stars that 
could be affected by these effects. 

In the calculation we use isochrones from the Padova group 
(Girardi et al. 2002), for LogAge(yr) values of 6.6, 6.7, 6.8, 6.9, 
7.0, 7.2, 7.4, and 8. The isochrones are first shifted in all CM di- 
agrams to account for the colour excess and distance moduli cal- 
culated for the probable unevolved MS members. Furthermore, 
the colour calibration used to translate the theoretical tracks to 
the observational plane present offsets with respect to the ZAMS 
of Schmidt-Kaler (1982) or Kenyon & Hartmann (1995), which 
coincide quite closely. When using the isochrones for age esti- 
mates, these offsets are previously corrected. 
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Fig. 2. Example of selection of probable unevolved MS members. The left panel shows the tentative selection of upper MS members 
in the V, (B - V) diagram of NGC 2367 (see text). In the right panel, the quantity V - 3.1 x E(B - V) is plotted for these stars versus 
My. The lines indicate published values of the distance modulus: Vogt & Moffat 1972 (continuous line) and Kharchenko 2005 
(broken line). Circles are stars selected from the plot as probable, unevolved MS members. Triangles and crosses are stars selected 
as post-MS members. Crosses mark the stars used to estimate the post-MS age of the cluster (see text) 



For each involved star, the difference in colour with respect 
to every isochrone is measured in the four CM diagrams avail- 
able. These differences are then either interpolated or extrapo- 
lated to the age value that would fit the star position exactly. 
This provides a nominal age value for each star, and all these 
values are then averaged to give the post-MS age estimate for the 
cluster. The average values and their rms deviations are listed in 
Table 3, together with other clusters parameters, and compared 
to published values. 

3.2. PMS membership 

We distinguish between stars for which a colour excess can be 
calculated in the CC diagram and stars that do not have this pos- 
sibility, either because of lacking U measurement, the most com- 
mon case, or because the shifting in the CC diagram does not 
provide any valid intersection with the reference lines. 

For stars with (U-B) measurement, the calculation proceeds 
as explained above. These stars were also compared with the 
PMS isochrones, taken from three different models (D'Antona 
& Mazzitelli 1997; Palla 1999; Siess et al. 2000, respectively re- 



ferred to in the following as D97, P99, and S00). The (U - B) vs 
(B - V) relation for these isochrones was the same as the one used 
for the ZAMS (Kenyon & Hartmann 1995). The colour excesses 
were therefore the same, whereas the distance moduli were dif- 
ferent, depending on the model and on the age of the isochrone. 
We used model isochrones of ages from 1 to lOMyr in steps 
of 1 Myr for the S00 models and 2Myr for D97 and P99 mod- 
els. Therefore, we obtained up to 20 possible DM values, from 
10 S00 isochrones, 5 P99 isochrones, and 5 D97 isochrones. A 
star was then assigned as a PMS member if, in addition to the 
colour excess requirement for PMS candidates (Sect. 3.1.2), one 
or more of these DM values agreed within the errors with the av- 
erage value of the probable MS unevolved members. Obviously, 
this membership assignment wss done simultaneously in the two 
CM diagrams V, (U - B) and V, (B - V). 

For stars without a valid (U-B) value, the CC diagram can- 
not be used, and we do not have calculated values of colour ex- 
cess and distance modulus to compare with reference lines. For 
these stars, only PMS membership was investigated. To assign 
PMS membership to these stars, we assumed every star was af- 
fected by the colour excess calculated for the probable unevolved 
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MS members. To compare stars with isochrones in the V, (V—R) 
and V, (V - I) CM diagrams, we used reddening slopes for ad- 
equate achieving of a good ZAMS fitting of the probable un- 
evolved MS members in the respective diagrams. These values 
for the reddening slopes are in good agreement with the values 
from a standard reddening law (Cardelli et al. 1999). The PMS 
isochrones are shifted with these values in the CM diagrams, and 
the comparison with the data gives a distance modulus for every 
star relative to each isochrone. A star has a membership assign- 
ment (with respect to an isochrone) if the obtained distance mod- 
ulus coincides, within the errors, with the average value for the 
probable MS unevolved members. In this procedure, a star can 
have a membership assignement in one, two, or all three CM 
diagrams involved. 

To clarify this procedure of member selection, we show 
examples in the CM diagrams of NGC3293 in Fig. 5 of stars 
with different qualities of membership assignment. The ellipses 
around each point reproduce the errors explained above, used for 
comparison: squared quadratic sums of the uncertainties plotted 
in Fig. 1, and rms deviation of the average colour excess and 
distance modulus for the probable unevolved MS members. The 
isochrones plotted in this example are from Siess et al. (2000). 



In both methods of membership assignement, the errors used 
in the comparisons are the squared quadratic sums of the uncer- 
tainties plotted in Fig. 1, and the rms deviation of the average 
values for the probable unevolved MS members. As can be ap- 
preciated in the example of Fig. 5, the errors are usually larger 
than the separation between two consecutive isochrones. As a re- 
sult, every star can be a candidate member with respect to several 
PMS isochrones. The average value of the corresponding ages is 
taken as the age of the PMS candidate. The average value of the 
ages for all assigned members and its rms value are adopted as 
the age and uncertainty of the PMS sequence of the cluster, as 
obtained from the particular CM diagram and model being com- 
pared. 

Two more points have to be covered before finally select- 
ing the candidate PMS members. First, the procedure described 
above for analysing the membership of stars without a valid 
(U-B) index was also applied to stars with a valid (U-B) value, 
and qualified as reddening candidates (see the last paragraph in 
Sect. 3.1.2). For stars without a (U-B) index, the average colour 
excess of the probable unevolved MS members was assumed, as 
explained above. For reddening candidates, their own calculated 
color excess value was used instead. In this way, PMS mem- 
bership assignment can be established for a star with respect to 



8 Delgado et al.: DAY-I Catalogue of PMS Stars 
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several combinations of three CM diagrams, and even with re- 
spect to four CM diagrams. This is indeed the case of a star with 
an eventual assignment based on its calculated colour excess and 
distance modulus (as mentioned above, this obviously means si- 
multaneous assignment in the V, (U - B) and V, (B - V) CM di- 
agrams), and also assigned as a member from its location in the 
V, (V - R) and V, (V - I) CM diagrams. 

Second, the procedure for membership assignment may lead, 
for some stars with a valid (U — B) value to an assignment both 



as MS and PMS members. This will happen for stars located in 
the region of the CM diagrams where the ZAMS and the oldest 
PMS isochrones merge. In these cases, we adopted the PMS as- 
signment for stars with spectral types later than AO (My =1.3 on 
the ZAMS) and the MS assignment for stars earlier than this. 

The final selection of candidate PMS members included only 
stars whith a membership assignment in at least three CM dia- 
grams. Furthermore, to achieve a conservative estimate, we con- 
sidered a further level of membership quality, with two more 
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Fig. 4. Example of the colour excess calculation in the (U - B), (B - V) diagram. Three stars in NGC 3590 with different types of 
solutions are plotted, together with the ZAMS (continuous line). Broken lines represent the shifts with the slope E(U-B)/E(B-V) = 
0.72. Ellipses represent the photometric uncertainties of the respective points, plotted in Figure 1. 



restrictions. A star is definitely considered as a possible PMS 
member if, in addition to being assigned in three CM diagrams, 
(i) its photometric errors (PSF errors) are below 0.05 in all five 
UBVRI bands and (ii) the ages obtained from all three (or even- 
tually four) CM diagrams coincide within a range of 1 Myr. The 
resulting assigned members are used to calculate the PMS age, 
as explained, and are the stars listed in our final catalogue. 



The results of the selection procedure are shown in the CM 
diagrams of the clusters. Figures 6 to 10 show the V, (B - V) 
CM diagrams for all clusters, compared to both post-MS (G02) 
and PMS (S00) isochrones. As mentioned before, three PMS 
model isochrones are used to calculate this membership assign- 
ment (D97, P99, and S00). Three different assignments and ages 
are therefore obtained for the PMS candidates in each clus- 
ter. In Figures 6 to 10, we only plot the candidate PMS mem- 
bers with respect to S00 isochrones. In our final catalogue, 
we include assigned members with respect to any of the three 
model isochrones. In this respect, we have to mention that D97 
isochrones cover only the faintest part of the CM diagrams, and 
they do not provide PMS membership assignment for stars ear- 
lier than GO. 



3.3. Astrometry and 2MASS infrared data 

The pixel coordinates of our catalogue were transformed 
to equatorial coordinates with the IRAF tasks ccmap and 
xy2rd, using matching files with stars in common between our 
fields and corresponding fields in the 2MASS data base. The 
matching provides equatorial right ascension and declination 
for all the stars observed, and JHK photometry for 15850 
stars, while 68% of them have photometric quality flag A 
to D in all three JHK bands (see the catalogue description at 
http : //www . ipac . caltech. edu/2mass/releases/all sky/doc/). 

The values from the 2MASS catalogue were transformed 
to the system of Bessel & Brett (1988) with relations given by 
Carpenter (2001). These values are also included in the final cat- 
alogue. 

3.4. Additional remarks 

Cluster membership assignments derived from kinematic studies 
have always been considered more reliable than those obtained 
from photometric analysis. However, the determination of prob- 
able cluster members from analysing their location in the photo- 
metric diagrams has an ample and successful history, ever since 
its origin around the middle of the past century. We have prof- 



10 



Delgado et al.: DAY-I Catalogue of PMS Stars 




0.5 1 1.5 0.5 1 0.5 1 1.5 

B-V V-R V-I 



Fig.5. Example of PMS members selection in the CM diagrams of NGC3293. Star # 54, non member; #1231, member in the 
V, (B - V) diagram; #334, member in three CM diagrams, with different ages; #2592, member in all three diagrams, with age 
differences among them below 1 Myr. The axes of the ellipses around each point reproduce the squared quadratic sums of final 
uncertainties (Fig. 1) and the rms deviation of cluster parameters for the probable unevolved MS members (Table 3). Isochrones of 
ages 1 to 10 Myr from Siess et al. (2000) are plotted 



ited from this experience, and with addition of new PMS evo- 
lutionary models, developed a methodology based on analysis 
of colour-colour and colour-magnitude diagrams to obtain a re- 
liable sample of PMS cluster candidate members. This method 
was designed and developed in previous works and applied to 
several clusters (Delgado et al. 1998, 2000, 2004, 2006). Here, 
we use it to study a cluster sample, with homogeneous U BVRcIc 
photometry, obtained at the same telescope and with the same 
instrumentation for all of them. The clusters are located mainly 
in the Carina complex of star formation (Alfaro et al. 1992), 
and they could therefore contain valuable information on the pri- 
mary evolutionary stage of stars, inside clusters that form part of 
a larger region of coherent star formation in our Galaxy. The 
study of clusters in a wide age range, from 4 to 30 Myr is of 
interest for analysing the possible sequential star formation and 
the resulting age spread inside the complex. Discussion of the 
age properties and of the members' spatial distribution and NIR 
properties is postponed to our next paper. 

Some different sources of uncertainty that, among others, 
might be affecting our membership results are: 

1. The (U - B), (B - V) relation for PMS stars, may differ 
from the one for MS stars, and the reddening law for PMS stars, 
where unquantified amounts of remaining circumstellar dust and 
gas may be present, can also be different from the reddening re- 
lations describing extinction for MS stars (Delgado et al. 1998). 

2. The values of the reddening slopes may indeed change be- 
tween clusters (Sagar & Cannon 1997; Guetter & Turner 1997). 



3. The eventual presence of high reddening makes it harder 
to distinguish PMS members and non-member stars (Chini & 
Wargau 1990). 

4. The uncertainties in the theoretical isochrones, which arise 
both from the model computation and from their translation to 
the observational plane as well. 

It is not possible to evaluate the total uncertainty in the fi- 
nal classification due to these factors. We do, however, have ar- 
guments to support its reliability. First, we recall the restrictive 
conditions for the stars to be selected as members, which re- 
quire not only a simultaneously consistent location in all CM di- 
agrams, but also a good photometric accuracy in all five UBVRI 
bands. 

On the other hand, we have performed some follow-up spec- 
troscopic studies of previously observed clusters, which provide 
results consistent with the photometric classification (Delgado et 
al. 1998, 1999, 2004). In the case of NGC2362 (Delgado et al. 
2006), a satisfactory agreement is reached between the described 
membership assignment and expectations from independently 
observed properties of the stars, such as Ha emission and X- 
ray activity. We think that the effects of all the above-mentioned 
sources of uncertainty is decreased in the presence of not very 
high reddening, as is decidedly the case for NGC 2362, and also 
- although not so favourable - as for the clusters in the present 
sample. 

Finally, we recall the comparison of our membership assign- 
ments with those by Baume et al. (2003) for NGC 3293, which 
shows a high degree of agreement: 89% of their 292 proposed 
members, also included in our observations, have some assign- 



Delgado et al.: DAY-I Catalogue of PMS Stars 



11 







Ill 
NGC 2367 


i 

1 


ill 
l Collinder 228 


o 

1— 1 


- 1 


If 






V 










LO 
H 


- 


v^H. ; : -J \ . 
■ \ % *?k 

/v-V/-: '.. \* ■ •\ 




3S6P a' \ r ^ 

' \ * \ - \ 


o 

oi 


— 




i 


\\\ 
1 1 \ \ 1 \ 






0.5 1.5 
B-V 




0.5 1.5 

B-V 



Fig. 6. V, (B - V) CM diagrams of NGC 2367 and Collinder 228, two of the youngest clusters in our sample. Circles represent the 
probable unevolved MS members. Crosses are the post:MS probable membrs selected to for the quantitative age estimate. Large dots 
represent PMS members in respect to any of the three model isochrones used, and selected from three CM diagrams. Squares are 
PMS members assigned in all four CM diagrams. Small dots are the remaining observed stars. The ZAMS line, post-MS isochrones 
by G02, for LogAge(yr) 6.6, 7.0, and 7.4, and S00 PMS isochrones for LogAge(yr) 6 and 7, are plotted for reference. 



ment by our procedure (either MS member or PMS member in 
at least one CMD). 

4. Catalogue and conclusions 

Multicolour UBVRcIc photometry has been obtained for 26962 
stars in the fields of 10 young clusters, 9 of them located in 
the Carina star formation complex of the Milky Way. Ages of 
the clusters, as estimated from their MS and post-MS probable 
members, range between LogAge(yr) 6.6 and 7.5. Matching of 
the pixel coordinates with the 2MASS catalogue provides abso- 
lute coordinates for DAY-I catalogue and JHK photometry for 
the stars in common. 

Analysis of the joint evidence from the CC and CM dia- 
grams, together with the comparison of the observational CM 
diagrams to theoretical PMS and post-MS isochrones, led to es- 
timates of both MS and PMS membership for individual stars 



in every cluster field. A total of 842 stars have been selected as 
PMS members in the different clusters, with respect to at least 
one of the three sets of isochrone models used to estimate PMS 
membership. For some stars, membership assignment and cor- 
responding age values were obtained according to these three 
models. 

For the first time, a catalogue of probable PMS cluster 
members has been elaborated. This catalogue is available at 
the CDS (Centre de Donnees astronomiques de Strasbourg) 
at http://cdsweb.u-strasbg.fr/CDS-f.html, and it con- 
tains the basic information from the calculations and analysis de- 
scribed here. The column contents is as follows. Column 1: a 15- 
character string with the star identification. The first 9 characters 
are the IAU cluster number, and the last five characters are the 
identification number in our photometric files. Columns 2 and 
3 list the right ascension and declination in degrees for Epoch 
2000; Columns 4-8 the optical photometric values V, (U — B), 
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(B - V), (V - R), and (V - I); Columns 9-11 the NIR photomet- 
ric values K, (J - H), and (H - K); Column 12 the membership 
indication: 1, PMS member assignments in three CM diagrams; 
2, PMS member assignments in all four CM diagrams. Column 
13 indicates absorption Ay = 3.1 x E(B - V) when available. 
Columns 14-16 show PMS ages with respect to the S00, P99, 
and D97 models, respectively. In Table 4 we list the first 10 lines 
of the catalogue as an example, which correspond to stars in 
NGC 2367. 

The catalogue is meant as a starting point for more de- 
tailed studies of PMS stars, either individually or as members of 
hierarchically-structured stellar systems. In particular, we con- 
sider it as a basis for follow-up spectroscopic analysis and even- 
tual observations in other wavelength ranges, mainly in the NIR 
range. 
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Fig. 10. The same as Figure 3, for the clusters NGC 4463 and Trumpler 18. 
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